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I .  introduction 

During  the  3-year  period  of  Contract  FI 9628-83-K-0037 ,  we  have 
carried  out  experimental  studies  of  ion-molecule  reactions  and 
electron-ion  recombination  processes  using  drift-tube  and  microwave 
afterglow  techniques. 

The  studies  of  ion-molecule  reactions  focussed  on  oxidation 
reactions  of  metal  ions  and  three-body  association  of  atmospheric 
diatomic  ions.  Specifically,  we  investigated: 

a)  Reactions  of  zirconium  ions  with  oxygen,  nitric  oxide,  and  carbon 
dioxide.  The  details  of  the  measurements  and  the  results  of  this  work 
are  described  in  an  article  published  in  The  Journal  of  Chemical 
Physics  (Appendix  I) 

b)  The  dependence  on  temperature  of  three-body  association  of 
diatomic  atmospheric  ions  with  atmospheric  gases  had  been 
investigated  by  us  in  the  previous  contract  period.  A  theoretical 
model  of  the  reaction  mechanisms  underlying  the  association  has  been 
developed  and  a  publication  describing  this  work  has  been  accepted  by 
The  Journal  of  Chemical  Physics  (Appendix  II). 

Our  experimental  studies  of  electron-ion  dissociative 
recombination  involved  both  simple  diatomic  ions  (N0+)  and 
clustered  species  (04+).  The  principal  emphasis  was  on 
determining  the  dependence  of  the  recombination  coefficients  on 
electron  temperatures  via  microwave  heating  of  electrons  in 
recombination  controlled  afterglows.  As  part  of  this  work  extensive 
computer  models  of  the  electron  energy  loss  and  gain  in  afterglow 
plasmas  containing  small  concentrations  of  molecular  additives 
(inelastic  losses)  have  been  constructed  and  have  been  applied  in 
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the  analysis  of  the  N0+  and  04+  recombination  data.  The 
present  measurements  and  analysis  for  N0+  recombination  have 
removed  the  discrepancy  with  trapped-ion  data  and  will  be  submitted  to 
Physical  Review  A  for  publicatici  shortly  (the  abstract  is  included 
with  this  report,  see  Appendix  II). 

II .  Research  published  during  the  grant  period 

The  texts  of  the  following  papers  appear  in  the  Appendix  of  this 
report: 

Ion-molecule  reactions: 

Rate  coefficients  for  the  oxidation  of  zirconium  ions  with 
oxygen,  nitric  oxide,  and  carbon  dioxide.  S.  Dheandhanoo,  B.K. 
Chatterjee,  and  R.  Johnsen,  J.  Chem.  Phys.  83,  3327  (1985) 

Temperature  dependence  of  association  of  diatomic  ions  in 
diatomic  gases,  R.  Johnsen,  J.  Chem.  Phys.  xx,  xxxx  (1986) 

III.  Research  results  in  preparation  for  publication: 

Electron  ion  recombination: 

Electron  temperature  dependence  of  the  recombination  of 
electrons  with  N0+  ions.  J.L.  Oulaney,  M.A.  Biondi,  and  R. 

Johnsen.  Phys.  Rev.  A,  to  be  submitted. 
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Appendix 


Text  of  titles  and  abstracts  of  papers  published,  accepted  for 
publication,  or  to  be  submitted  shortly  in  scientific  journals  during 
the  period  of  this  contract  are  given  on  the  following  pages. 


4 


i 


Rate  coefficients  for  the  oxidation  reactions  of  zirconium  ions  with  oxygen, 
nitric  oxide,  and  carbon  dioxide 

S.  Dheandhanoo,  B.  K.  Chatterjee,  and  R.  Johnsen 

Deportment  of  Physics  and  Astronomy.  University  of  Pittsburgh.  Pittsburgh,  Pennsylvania  15260 
(Received  17  May  1985;  accepted  20  June  1985] 

Rate  coefficients  for  the  oxidation  reactions  of  zirconium  ions  with  oxygen,  nitric  c.tide,  and 
carbon  dioxide  have  been  determined  in  the  near-thermal  energy  range  (0.04  to  about  0.4  eV) 
using  a  selected-ion  drift  apparatus.  The  reactions  are  found  to  proceed  with  large  rate  coefficients 
[~  5|  —  10)cm  Vs  at  300  K).  At  elevated  ion  energies  the  reactions  with  NO  and  CO-  exhibit 
smaller  rates.  The  secondary  oxidation  reactions  of  ZrO  *  ions  with  the  same  reactants  have  rate 
coefficients  about  two  orders  of  magnitude  smaller  than  those  for  the  first  oxidation  steps.  The 
discharge  ion  source  used  to  produce  the  zirconium  ions  in  these  experiments  is  described  in  some 
detail. 


I.  INTRODUCTION 

While  gas-phase  oxidation  reactions  of  metal  ions  with 
atmospheric  gases  are  of  considerable  applied  interest,  com¬ 
paratively  little  experimental  results  are  available  In  an  ear¬ 
lier  series  of  measurements'^  we  examined  reactions  of  the 


type 

M*  +  O,— MO  *  +  O 

111 

and 

MO  “  4-  02—M02‘  +  O 

(2) 

for  ions  M  *  =  U  * ,  Th  * ,  and  Ti  ‘  and  found  that  the  first 
oxidation  step  generally  occurred  with  large  rate  coeffi¬ 
cients,  while  secondary  oxidation  was  usually  a  slow  process. 
We  have  extended  this  work  to  oxidation  reactions  of  zirco¬ 
nium  ions  with  02,  CO-,  and  NO  and  present  he;e  the  mea¬ 
sured  rate  coefficients  and  their  dependence  on  ion  energy  in 
the  near-thermal  energy  range. 

The  experimental  difficulty  in  such  experiments  lies  pri¬ 
marily  with  producing  the  ions,  a  problem  even  more  severe 
in  the  case  of  the  refractory  metals.  The  ion  source  developed 
for  this  work  may  be  applicable  to  a  variety  of  metals  and  is 
described  in  some  detail 

II.  EXPERIMENTAL  METHODS 

The  measurements  were  carried  out  using  a  selected  ion 
drift  apparatus  (SI DA),  which  except  for  the  Zr  ’  ion  source 
had  the  same  configuration  as  described  earlier. '  The  “addi¬ 
tional  residence  time”  method4  employing  programmed 
drift  fields  to  vary  the  ions’  residence  time  in  the  drift/reac¬ 
tion  region  was  used  to  measure  the  decay  of  Zr  *  tons  in  the 
presence  of  the  reactants,  from  which  the  rate  coefficients 
were  derived. 

It  proved  difficult  to  obtain  adequate  Zr*  ion  currents 
from  ion  sources  based  on  the  principle  of  evaporating  the 
metal  from  heater  substrates  with  subsequent  electron  im¬ 
pact  ionization,  a  method  employed  in  earlier  work. In  the 
case  of  zirconium,  such  sources  were  plagued  by  rapid  de¬ 
struction  of  the  heaters  (W  wire  or  indirectly  heated  tubes  of 
alumina  or  zirconia  were  tried  I  and  required  too  frequent 
replenishment  of  the  Zr  metal  samples.  A  more  satisfactory 
ion  source  was  built  in  the  form  of  a  magnetically  confined 
(Penning)  discharge,  shown  schematically  in  Fig.  1,  which 


was  operated  in  a  mixture  of  neon  and  zirconium  letrabro- 
mide  vapor.  It  consisted  of  a  cylindrical  discharge  chamber 
(5  cm  long.  2.5  cm  in  diameter),  formed  by  a  cylindrical, 
stainless  steel  anode  and  two  cathode  end  plates,  one  con¬ 
taining  the  ion  extraction  orifice  (0.25  cm  in  diameter5.  Mica 
gaskets  served  to  insulate  and  seal  the  anode  and  cathodes 
with  respect  to  each  other.  A  permanent  magnet  surround 
ing  the  anode  provided  the  axial  magnetic  field 

Approximately  1  g  of  ZrBr,  powder  was  placed  into  the 
discharge  chamber  and  evaporated  during  operation  of  the 
source  by  heating  the  entire  source  assembly  to  about  160  ’C 
To  obtain  a  stable  discharge,  neon  gas  was  passed  through 
the  source  at  pressures  of  0. 1  Torr  and  the  discharge  current 
was  adjusted  to  typically  0.5  to  1.0  mA.  After  about  100  h  of 
operation  the  ZrBr4  supply  was  exhausted  and  had  to  be 
replenished. 

In  addition  to  Zr*  ions  the  source  produced  large  ion 
currents  of  Ne*  and  Br*.  which  were  prevented  from  enter¬ 
ing  the  drift  region  of  the  apparatus  by  the  ion-selecting  mass 


«°»  inler  heater 


FIG  1  Schematic  diagram  of  the  zirconium  ion  discharge  soukc 


J  Chem  Phys  #3(7),  i  October  1985  0021  -9606/85/  193327-03S02  10  c  ’985  American  Institute  0*  Physics  3 327 

Th#  Vj.S.  Government  ts  authorized  to  reproduce  and  sen  thi*  report 
^•rmJtslon  for  further  reproduction  by  other*  must  be  obtained  from 
the  copyrleht  owner 


6 


3328 


Dheandhanoo.  Chatter  j®@,  and  Johnson  Oxidation  reactions  of  zirconium  ions 


filter  For  the  measurements  pulsed  operation  of  the  source 
was  desired.  Since  it  was  found  to  be  impractical  to  achieve 
pulsing  of  the  discharge  by  switching  the  high-voltage  sup¬ 
ply.  the  source  was  dc  biased  such  as  to  cut  off  the  Zr*  ion 
current  entering  the  drift  tube.  Application  of  a  short  (10/us 
longi  200  V  pulse  to  the  bias  network  then  produced  pulsed 
injection  of  ions  into  the  drift  tube. 

III.  EXPERIMENTS  AND  RESULTS 

Zirconium  ions  from  the  ion  source  were  injected  into 
the  drift  tube  containing  typically  0.5  Torr  of  helium  buffer 
with  small  admixtures  of  the  reactants  O,.  NO,  or  C02 
tabout  10  4  Tom.  The  ions  transmitted  through  the  drift 
tube  w  ere  found  to  consist  largely  of  the  Zr  *  parent  ions  and 
the  ZrO  *  products  produced  by  the  reactions: 


Zr* 

-  O,— ZrO*  +  0, 

(3) 

Zr* 

*  NO— ZrO*  +  N, 

<4| 

Zr‘ 

-  CO,— ZrO*  +  CO 

(5) 

As  several  minor  impurity  ions  were  present  in  the  mass 
range  of  interest,  the  isotopic  abundances  of  both  parent  and 
product  ions  were  examined  carefully  to  prevent  incorrect 
identification  of  the  ion  species.  At  higher  reactant  concen¬ 
trations.  further  oxidation  of  ZrO*  ions  led  to  the  appear¬ 
ance  of  ZrO,*  ions  in  the  mass  spectrum. 

After  having  established  that  the  oxidation  reactions 
constituted  the  dominant  loss  mechanism  of  Zr*  ions,  the 
rate  of  loss  of  Zr  *  ions  was  determined  using  the  additional 
residence  time  techniques  at  several  values  of  the  field-to- 
density  ratios  £  /.V,  which  controls  the  Zr  *  ions’  mean  ener¬ 
gy  The  well  known  expression  given  by  Wannier'  was  used 
to  infer  the  mean  ion  energy  from  the  measured  Zr*  drift 
velocities  The  energy  scale  shown  in  the  Figs.  2  through  4 
corresponds  to  this  energy  converted  to  the  Zr* /reactant 
center-of-mass  system  of  reference. 

The  measured  rate  coefficients  for  reactions  (3),  (4),  and 
|5l  are  shown  in  Figs.  2  through  4  as  function  of  the  Zr* 
center-of-mass  energy.  The  thermal  energy  (0.04  eV)  rate 


0  0.4  0.8  1.2 


<*cm>  <«V) 

FIG  2  Measured  rite  coefficients  for  the  reaction  of  zirconium  ions  with 
osvgen  as  a  function  of  the  mean  ion  energy  Data  point  at  004  eV  repre¬ 
sents  an  average  of  many  measurements 


<*cm> 

FIG.  3.  Measured  rate  coefficients  for  the  reaction  of  zirconium  ions  with 
nitric  oxide  as  a  function  of  the  mean  ion  energy  Data  point  at  0  04  eV 
represents  an  average  of  many  measurements 

coefficients  are  close  to  the  limiting  Langevin  rate  coeffi¬ 
cients  in  all  three  cases,  but  only  the  reaction  with  oxygen 
exhibits  an  energy  independent  coefficient  over  the  energy 
range  covered  by  the  experiments.  Fast  ion  molecule  reac¬ 
tions  exhibiting  rate  coefficients  near  the  Langevin  limit 
usually  show  little  or  no  dependence  on  ion  energy.  The  reac¬ 
tions  of  zirconium  ions  with  NO  and  C02  appear  to  be  an 
exception  to  this  rule,  which  is  somewhat  surprising.  It  is 
possible  that  these  reactions  proceed  through  an  intermedi¬ 
ate  reaction  complex,  whose  lifetime  decreases  with  increas¬ 
ing  ion  energy  leading  to  a  reduced  rate  coefficient  at  higher 
energies. 

The  binding  energy  of  ZrO  *  apparently  is  not  known 
from  other  measurements,  such  that  it  is  difficult  to  deter¬ 
mine  the  exoergicitics  for  reactions  (3),  (4),  and  (5).  The  oc¬ 
currence  of  reaction  (4)  indicates  that  the  binding  energy  of 
ZrO*  is  greater  than  that  of  nitric  oxide  (6.54  eV).ft  The 
isoelectronic  molecule  yttrium  oxide  has  a  binding  energy  of 
7.02  eV.6  Thus  a  similarly  large  binding  energy  for  ZrO* 
seems  plausible. 

Secondary  oxidation  reactions  of  ZrO*  ions  with  02, 
NO,  and  CO,  were  observed  to  occur  but  with  rate  coeffi¬ 
cients  two  to  three  orders  of  magnitude  smaller  than  those 


<  ecm  > 


FIG  4  Measured  rale  coefficients  for  the  reaction  of  zirconium  ions  with 
carbon  dioxide  aa  a  function  of  the  mean  ion  energy  Data  point  at  0  04  eV 
represents  an  average  of  many  measurements 
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for  the  first  oxidation  steps.  Since  it  was  found  difficult  to 
obtain  experimental  conditions  suitable  for  accurate  deter¬ 
minations  of  the  secondary  rate  coefficients  only  approxi¬ 
mate  rate  coefficients  were  derived  from  the  observed  abun¬ 
dance  of  the  secondary  product  Zr02*  in  the  ion  mass 
spectrum. 

For  the  reaction 

ZrO  '  +  O,— Zr02*  +O  |6| 

a  thermal-energy  (300  K)  rate  coefficient  of  (5  *2,|x  10 
cm'/s  would  be  required  to  explain  the  observed  ZrO;* 
|  abundance.  Corresponding  estimates  for  the  reactions 

ZrO*  ■+■  NO— ZrO,*  +  N  |7| 

and 

ZrO  -CO,— Zr02*  +  CO  |8) 

yielded  values  of  k  |7|  =  (2.5  *  ,a/s | x  10  12  cm'/s  and 
£I8|  =  ||  *,),txl0  12  cm  Vs.  While  the  energy  dependence 
of  these  reactions  was  not  studied  systematically,  the  few 
measurements  that  were  made  indicated  that  the  rate  coeffi¬ 
cients  are  likely  to  decrease  with  increasing  ion  energy. 

IV.  SUMMARY  AND  CONCLUSIONS 

The  measurements  confirm  the  expectation  that  oxida¬ 
tion  reactions  of  zirconium  ions  with  O,,  NO,  and  CO, 


should  be  rapid  processes  with  rate  coefficients  ciose  to  the 
limiting  Langevin  value.  The  energy  dependence  of  the  reac¬ 
tions  with  NO  and  CO,  appears  to  indicate  that  these  reac¬ 
tions  proceed  through  an  intermediate  reaction  complex 
rather  than  by  a  direct  mechanism. 

Secondary  oxidation  reactions  of  ZrO*  proceed  with 
very  small  rate  coefficients  and  hence  may  not  be  important 
in  applications. 
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Temperature  dependence  of  association  of  diatomic  ions  in  diatomic  gases 

Rainer  Johnsen 

Department  of  Physics  and  Astronomy.  University  of  Pittsburgh,  Pittsburgh,  Pennsylvania  15260 
(Received  18  April  1986;  accepted  24  June  1986) 

The  temperature  dependence  of  experimental  rate  coefficients  for  association  of  diatomic  ions 
in  diatomic  gases  is  compared  to  model  calculations  in  which  the  stabilization  of  the  initially 
formed  unstable  association  complex  is  assumed  to  occur  by  formation  of  a  “super  complex." 

It  is  shown  that  acceptable  agreement  between  the  model  and  experimental  data  on  four 
sample  reactions  can  be  obtained  with  a  minimum  of  plausible  assumptions.  The  agreement 
indicates  that  earlier  theoretical  models  may  provide  an  accurate  description  of  the  formation 
and  destruction  of  the  association  complex,  but  failed  to  reproduce  the  experimental  data, 
because  the  temperature  dependence  of  the  stabilization  step  was  inadequately  treated. 


I.  INTRODUCTION 

The  association  of  ions  in  gases,  that  is  a  process  of  the 

type 

A*  +  B  +  M— AB  •  +  M  (1 ) 

has  been  the  subject  of  many  experimental  and  theoretical 
investigations,  but  the  problem  of  reconciling  theoretical 
and  experimental  results  still  persists,  especially  in  regard  to 
the  dependence  of  the  association  rate  coefficients  on  tem¬ 
perature.  In  one  of  his  many  contributions  to  the  theory  of 
ion  association,  Bates'  referred  to  the  problem  as  an  "enig¬ 
ma,"  since  it  is  difficult  to  identify  faults  in  either  the  statisti¬ 
cal  theories  of  the  association  mechanism  or  the  experimen¬ 
tal  data 

Most  treatments  of  the  problem1'*'  have  focused  on  the 
first  step  in  the  association  mechanism,  formation  and  decay 
of  the  intermediate  unstable  complex.  The  second,  stabiliza¬ 
tion  step  is  either  assumed  to  proceed  with  a  temperature 
independent  rate  coefficient  or  an  essentially  ad  hoc  tem¬ 
perature  dependence  is  assigned  to  the  stabilization  step  in 
order  to  obtain  agreement  of  theory  with  experimental  data, 
a  clearly  unsatisfactory  procedure.  In  this  paper  it  is  suggest¬ 
ed  that  the  stabilization  of  the  reaction  complex  proceeds 
through  a  second  intermediate  complex,  called  the  super 
complex  It  is  shown  that  such  a  model  leads  naturally  to  a 
strongly  temperature  dependent  complex  stabilization  and 
that  is  capable  of  reproducing  experimental  data  for  diatom¬ 
ic  ion  association  in  diatomic  gases  with  acceptable  accura¬ 
cy.  The  basic  steps  in  the  association  mechanism  are  re¬ 
viewed  in  order  to  clarify  terms  and  to  state  the  assumptions 
on  which  it  is  based: 

The  rate  at  which  the  ions  A  ‘  are  converted  to  the  asso¬ 
ciated  ions  AB*  is  thought  to  be  limited  by  two  distinct 
mechanisms,  formation  of  an  intermediate  unstable  complex 
and  subsequent  collisional  stabilization.  If  either  or  both  A 
and  B  are  molecular  species  (we  exclude  the  special  case 
where  both  are  atomic)  a  fraction  of  the  energy  liberated  in  a 
collision  of  A  *  with  B  may  be  temporarily  absorbed  by  in¬ 
ternal  modes  of  the  AB  *  collision  complex,  i.e  . 

A  '  -  B  -(AB*  )*».  (2) 

where  the  symbol  *•  indicates  that  the  complex  has  suffi¬ 
cient  energy  to  redissociate  into  the  initial  particles: 


The  internal  dynamics  of  the  highly  excited  complex  is  ex¬ 
pected  to  be  quite  different  from  that  of  the  corresponding 
ground  state  molecule  and  it  is  not  obvious  that  a  definite 
geometric  structure  should  be  assigned  to  the  complex  Rap¬ 
id  exchange  of  energy  among  internal  rotations  and  vibra¬ 
tions  will  lead  to  highly  irregular  motions  in  the  potential 
well,  which  one  cannot  hope  to  describe  in  detail  Whether 
or  not  the  assumption  of  complete  randomness  t  "chaos"  )  is 
valid  for  a  particular  complex,  is  a  subject  of  much  current 
debate,  but  we  are  here  only  concerned  with  the  statistical 
properties  of  ensembles  of  such  complexes,  not  the  dynamics 
of  a  particular  complex.  If  it  is  assumed  that  both  A  *  and  B 
are  members  of  thermally  equilibrated  ensembles  of  particles 
at  a  common  temperature  T  and  that  the  unstable  reaction 
complex  has  no  modes  of  decay  in  addition  to  redissociation 
(e  g.,  ion-molecule  reactions  or  radiative  stabilization),  re¬ 
actions  (2)  and  (?)  will  lead  to  equilibrium  populations 

[  ( AB*  )**)/(  A  *  ]  =  [  B ) A: ( 2 ) / r .  (4) 

where  the  square  brackets  denote  concent  rations.  k\  2 )  is  the 
bimolecular  rate  coefficient  of  reaction  (2),  and  T  the  decay 
rate  of  an  ensemble  of  unstable  complexes  [  reaction  ( ? )  ]  in 
thermal  equilibrium  The  definition  of  T  requires  some  ex¬ 
planation:  If,  in  a  hypothetical  experiment,  after  the  equilib¬ 
rium  between  reactions  ( 2 )  and  (  ? )  have  been  attained,  re¬ 
action  (2)  is  made  to  stop  by  some  means  at  time  the 
density  nit)  of  unstable  complexes  will  initially  decay  as 

------  =  —  T«(/)  at  /=:/„.  (5) 

dt 

which  defines  T  at  t  =  r,,  only  ,  but  T  itself  may  depend  on 
time,  since,  for  instance,  complexes  with  large  total  energy 
may  decay  first  Thus,  the  solution  of  Eq.  ( 5)  is  not  necessar¬ 
ily  an  exponential  function  and  it  is  not  permissible  to  identi¬ 
fy  the  reciprocal  of  T. 

r=  r  1  ;  6) 

with  the  average  lifetime  of  the  ensemble  of  complexes  The 
properly  defined  average  lifetime  is  given  by 


and  may  differ  greatly  from  T  1  or  may  even  be  infinite 
This  formal  point  has  been  made  earlier  ‘  There  is  nothing 
wrong  with  identifying  T  as  a  reciprocal  lifetime  in  qualita- 
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live  discussions,  but  erroneous  conclusions  may  be  drawn, 
for  instance,  if  one  inserts  r  as  defined  by  Fq  ( 6 )  into  equa 
lions  or  inequalities  w  hich  compare  it  for  instance,  to  other 
characteristic  times,  such  as  the  poorly  defined  "randomiza¬ 
tion  time"  of  a  particular  complex.  The  time  available  for 
accessing  specific  states  of  a  particular  complex  may  be  ter¬ 
minated  by  collisions  (see  below)  that  eliminate  it  from  the 
ensemble  but  not  by  the  faster  decay  of  other  complexes  The 
confusion  arises  only  if  one  fails  to  distinguish  carefully  the 
fate  of  particular  complexes  from  that  of  the  ensemble  to 
which  they  belong  In  general,  the  equilibrium  of  reactions 
(  2  i  and  t 3 )  w  ill  be  slightly  perturbed  by  stabilizing,  energy 
transfer  collisions  of  the  complex  with  the  ambient  gas  M, 

(AH'  '*•  a  M  »(AB-  )*  q  M.  (8) 

where  ihe  single  asterisk  is  meant  to  indicate  that  the  AB  * 
ion  may  still  be  in  an  excited  state  but  has  insufficient  energy 
to  dissociate  Reaction  (8)  is  meant  to  describe  the  net  effect 
of  energy  removing  collisions  including  possible  colltsional 
reexcitation  of  highly  excited  but  stable  ions  to  states  above 
the  dissociation  limit  It  does  not  necessarily  represent  a  sin¬ 
gle  atomic  collision  process  and  no  elementary  cross  section 
should  be  assigned  to  it  The  experimental  data  on  the  den¬ 
sity  dependence  of  three-body  association  indicate  that  if  is 
valid  to  represent  the  stabilization  by  an  effective  bimolecu- 
lar  rate  coefficient  in  the  low-pressuie  limit,  but  they  yield 
little  information  concerning  its  magnitude  or  dependence 
on  temperature. 

The  low  pressure  limit  of  three-body  association  is  de¬ 
fined  by  the  condition 

k(  8 )  [  M  |  <T,  (9) 

which  isequivalent  to  the  assumption  that  the  equilibrium  of 
reactions  ‘  2  '  and  (3)  is  only  negligibly  perturbed  by  stabi¬ 
lizing  collisions.  In  the  low-pressure  limit  the  overall  rate 
coefficient  for  three-body  association  k  may  then  be  written 


as 

A  t 

i  /  ;A( 8). 

(10) 

w  here 

1  ’  /  , 

A  i  2  )  /  r 

(11) 

is  (he  tempo ature-dependent  equilibrium  constant  for  reac- 
lions  i  2  '  and  (  3 )  Inequality  ( u )  also  implies  that  the  ener¬ 
gy  distribution  of  the  unstable  complexes  is  only  negligibly 
perturbed  by  (he  stabilizing  collisions  At  higher  densities  of 
M.  bul  lit  I  in  the  low  -pressure  limit,  (his  approximation  may 
•vcome  invalid  and  A (  8 )  will  no  longer  be  independent  of 
ga-  density .  i  e  .  it  is  not  a  proper  rate  coefficient  It  is  for  this 
cavil  ilia!  a  is  not  entirely  correct  to  extend  Fq  (10)  to 
higher  gas  densities  m  the  form 

s  I  /  -A  I  x>r.'(  I  -  A(8i|  \f  I  ),  (12) 

a-  is  s.  inclimes  done 

In  olds’  to  compare  the  low-pressure  ute  coefficient 
cak  u'ateil  liom  I  q  :  10)  with  experimental  data,  it  is  neces¬ 
sary  -  bi.ur,  the  factors  f  (  /)  and  A  (  8  )  Several  authors 
'i‘Vi  ,  ik  mated  I '(  /  *  by  methods  based  on  phase  space  the- 
’i  v.  isKKd  'heory.oi  extensions  of  this  approach  1  v  Of 
parti,  alai  mieiest  is  the  temperature  dependence  of  H  I  ), 
for  which  the  following  useful  "rule  of  thumb"  has  been 


derived  (see  Ref  3,  bq  43)  A '(  7  )  depends  on  7' in  the  form 

»T  7  )  l  <  r„)<  /  //;,)  *.  (13) 

where 

x  1R  72.  i  14) 

and  R  represent  the  rotational  degrees  of  freedom  m  ihe 
reactants  A  '  and  B  It  A  and  B  are  diatomu  »  2  for 

instance  Vibrational  degrees  ot  freedom  may  be  ncglec  led  so 
long  as  the  vibrational  sp.n  mg  in  ihe  reactants  is  small  com 
pared  to  thermal  energy  I  he  simple  i  ule  has  a  simple  inter 
pretation  The  decay  rate  I  is  proportional  to  the  produv  t  of 
the  internal  partition  functions  ot  its  decay  products  Kach 
rotational  degree  of  freedom  .ontr\butes  a  fac  tor  /  to  the 
dr  pendence  of  T  on  / 

If  one  assumes  that  its  rale  'im.iiii  tor  stabilization 
A <  b )  is  independent  on  /  u  I  >llows  ,u  once  that  A  should 
depend  on  7  as 

A  A  (  /„)(  /  /  i  1  I  |X  i 

in  ihe  case  of  dialomii  re.ii lanis  i  2  As  will  be  seen  lalcr. 
this  simple  prediction  is  supported  by  experimental  data 
only  for  some  but  not  all  reactions  One  often  finds  that  the 
observed  dependence  on  /  is  stronger  than  that  given  by  Eq 
I  I  X  ).  and  it  has  been  common  practice  to  replace  a  by  a  •  f> 
in  that  equation  The  additional  7  *  dependence  is  some 

limes  ascribed  to  the  variation  of  k  ( 8 )  with  7  I  he  practice  is 
not  only  an  unsatisfai  lory ,  ad  hoi  addition  to  theory .  but  it  is 
also  clear  that  the  functional  form  of  the  T  dependence  of 
A  ( 8 )  cannot  be  correct  at  all  temperatures  The  Langcvin 
collision  rate  constant  should  be  an  upper  bound  for  A(  8 ), 
which  would  be  exceeded  if  the  T  A  dependence  remained 
valid  at  low  temperatures 

II.  COMPARISON  OF  EXPERIMENTAL  DATA  AND 
THEORY 


The  set  of  reactions  chosen  for  comparing  theory  and 
experiment  involves  diatomic  ions  and  diatomic  neutrals. 
The  specific  reactions  are: 


NT 

a  N.  q 

N. 

•N/  - 

NT, 

(  16) 

OT 

a  O.  a 

O, 

-( )4  * 

(),, 

(17) 

( ).  ‘ 

i  N,  a 

N , 

.O  T  N 

•  •  N„ 

(  18) 

NO 

'  *-  N-  • 

.  V 

.NO- 

V  ;  N 

1  1 9 ; 

All  four  react 

lions 

,  involve 

reactants  of  similar 

molecular 

onsiaiil' 

>.  but  th 

e  si; 

ibihly  ot 

the  clustered  ions 

decreases 

from  about  1  eV  for  the  first  reaction  to  about  0.2  eV  for  the 
last  one  Much  of  the  relevant  information  has  been  com¬ 
piled  m  Table  I  which  also  contains  references  to  the  pub¬ 
lished  data  and  power-law  fits  to  the  experimental  results 
Reaction  ( lb)  has  been  measured  by  many  authors  in  addi¬ 
tion  to  the  referenced  data  Ihe  results  are  in  acceptable 
agreement  (on  the  order  of  30C£  deviations  over  the  tem¬ 
perature  ranges  covered)  and  the  temperature  exponents 
range  from  1.7  to  2  14 

The  rate  coefficient  of  reaction  (17)  was  once  believed" 
to  exhibit  a  peculiar  maximum  ai  low  temperatures  Later 
experimental  data"’  did  not  show  this  maximum  and  Fergu¬ 
son21  has  pointed  out  that  such  a  behavior  should  never  be 
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TABLE  I  For  reactions  ( lf>)-(  19)  (numberedastn  Sec  II).  the  table  lists 
experimental  values  for  the  threebod>  rate  coefficients  ay  300  K.  the  expo¬ 
nents  i  assuming  a  pcmer-lavv  temperature  dependence  4(7*)-  l  7*  / 
300 1  ind  references  to  the  espenmental  data  The  last  two  columms 
gisethemeasu  red  bond  energies  meV  of  the  product  cluster  ions,  and  refer¬ 
ences 


Reaction 

300  K  )  (cmVs) 

X 

Ref 

ieV) 

Ref 

5  . 

10 

2  u 

10 

1  0 

14 

r 

4  • 

10 

2  43 

11 

0  42 

1? 

IS 

1  ■ 

10  " 

3.2 

10 

0  22 

16 

l» 

}  - 

10  “ 

4  4 

10 

0  18 

13 

observed  There  are  few  measurements  of  the  remaining  two 
reactions,  but  one  might  add  that  a  strong  temperature  de¬ 
pendence  of  reaction  (19)  U~5)  has  also  been  inferred 
from  ionospheric  data  * 

As  is  clear  from  Table  I,  except  in  one  case  the  observed 
values  of  x  differ  from  t  -  2.  the  value  expected  from  the 
above  mentioned  rule  of  thumb.  Recently,  two  attempts 
have  been  made  to  rationalize  this  departure  from  simple 
expectation  Bates,'  in  an  attempt  to  explain  the  T  depen¬ 
dence  of  reaction  ( 17).  made  the  tentative  suggestion  that 
the  internal  states  of  the  04  *  unstable  complex  might  not  be 
fully  “randomized"  at  the  higher  temperatures.  Although 
his  estimates  of  the  additional  T  dependence  due  to  this  ef¬ 
fect  appeared  plausible,  he  later*  withdrew  his  suggestions  as 
untenable  in  the  low  -pressure  limit  of  the  reaction.  His  later 
position  is  correct.  As  discussed  in  Sec.  I.  departures  from  a 
thermal  equilibrium  distributions  in  the  complex  cannot  be 
invoked  in  the  low-pressure  limit.  In  his  later  paper,*’  Bates 
also  suggested  a  different  reason  why  some  vibrational  states 
in  the  complex  may  not  be  excited  at  all  temperatures,  but 
without  detailed  further  calculations  it  is  difficult  to  estimate 
the  magnitude  of  the  effect  In  one  of  his  earlier  papers 
Bates'  calculated  partition  functions  for  the  04'  complex 
(see  Table  I  in  his  paper)  assuming  either  absence  or  pres¬ 
ence  of  vibrational  excitation  of  the  constituent  molecules 
The  difference  between  the  two  cases  was  found  to  be  only 
about  25%,  which  makes  it  somewhat  unlikely  that  much 
additional  temperature  dependence  can  arise  from  vibra¬ 
tional  excitation 

Patrick  and  Golden*  presented  theoretical  calculations 
for  reactions  ( 16),  ( 17),  and  ( 18)  which  predicted  a  T  ~~ 
dependence  for  all  three  reactions.  In  the  case  of  reaction 
(18)  those  authors  compared  their  results  only  to  a  subset  of 
the  experimental  data  of  Dheandhanoo  and  Johnsen1" 
quoted  by  Speller  et  al. 2,1  Even  over  the  limited  range  of  tern  - 
peratures  where  the  comparison  was  made,  the  deviation  of 
the  data  from  a  T  2  dependence  is  quite  obvious.  Since  Pa¬ 
trick  and  Golden’s  results  for  reaction  ( 16)  gave  rate  con¬ 
stants  considerably  higher  than  the  experimental  data  by 
r  Bohringer  and  Arnold"  the  authors  suggested  that  the  ex¬ 

perimental  data  might  not  have  been  obtained  in  the  low- 
pressure  limn  Using  the  expression  in  Eq  ( 12 ).  they  show 
that  their  theoretical  results  can  be  made  to  agree  with  the 
experimental  data,  if  the  gas  pressure  in  the  experiment  was 
on  the  order  of  3  Torrat  T~65  K  Although  the  experimen¬ 
talists  did  not  clearly  state  the  gas  pressures  used  it  is  most 


unlikely  that  they  employed  such  high  pressure  at  low  tem¬ 
peratures:  From  the  figures  given  in  their  paper."  one  can 
reconstruct  the  conditions  of  the  experiment  and  finds  that 
the  pressures  used  at  the  lowest  temperatures  must  have 
been  less  than  0  I  Torr  To  give  a  simpler  argument:  If  the 
experimentalists  had  employed  a  gas  pressure  of  3  Torr,  the 
reaction  would  have  gone  to  completion  in  a  time  on  the 
order  of  several  nanoseconds,  w  hich  is  orders  of  magnitude 
shorter  than  practical  for  the  experimental  method  used. 
Thus.  Patrick  and  Golden's  reinterpretation  of  the  experi¬ 
mental  data  is  not  conv  incing. 

In  the  following  section,  the  disagreement  between  ex¬ 
periment  and  theory  is  seen  as  resulting  from  the  tempera¬ 
ture  dependence  of  the  complex  stabilization  rather  than 
from  deficiencies  in  the  calculations  of  l'(  T ). 

III.  STABILIZATION  VIA  FORMATION  OF  A  SUPER 
COMPLEX 

Since  the  temperature  dependence  of  I’(  T)  [  Eqs.  i  1 3 ) 
and  (14)]  seems  to  be  well  established  theoretically .  the 
source  of  the  discrepancies  between  theory  and  experiment 
may  be  found  in  the  temperature  dependence  of  the  stabiliza¬ 
tion  step  [  Eq  ( 8 )  ] .  The  detailed  mechanism  of  energy  re¬ 
moval  from  the  initially  formed  unstable  complexes  is  likely 
to  be  complicated  and  possibly  involves  several  collisions 
with  background  gas  molecules.  One  may  hope  that  the  net 
effect  of  the  energy  exchanging  collisions  can  be  described 
phenomenologically  by  assuming  the  existence  of  super 
complexes  consisting  of  the  unstable  complex  and  an  addi¬ 
tional  third-body  molecule.  The  idea  is  not  new  :  Herbst* 
applied  it  to  several  reactions  in  which  helium  was  the  third 
body  and  showed  that  this  model  leads  to  a  temperature 
dependent  stabilization  efficiency.  Formation  of  super  com¬ 
plexes  with  molecular  third  bodies,  the  case  of  interest  here, 
should  be  far  more  likely  than  with  weakly  interacting  rare 
gas  atoms,  such  as  helium.  After  all.  the  existence  ot  stable, 
higher-order  aggregate  ions  with  nitrogen  is  known  from 
experiments.*"  but  one  rarely  observes  ions  clustered  to  heli¬ 
um.  The  super  complexes  are  assumed  to  be  formed  in  the 
reaction 

( AB  ’)••-*-  M  —( ABM  '  )**  l  20! 

and  to  possess  two  modes  of  decay,  redissociation  into  the 
unstable  first  complex  and  M 

(ABM*  )**  •(  AB  )•*  -  M  ■  21  > 

or  dissociation  into  a  stable  (  AB '  )*  molecule  and  M 

(ABM*  )••--♦(  AB  '  )*  —  M.  (22) 

If  we  denote  the  bimolecular  rate  coefficient  of  Eq  i  20  by 
kt  and  the  unimolecular  rates  of  Eqs  (21)  and  i  22  >  by  T 
and  T  . .  respectively  ,  the  overall  rate  of  stabilization  k  ol 
the  first  complexes  is  given  by 

k.  —  k,  (r .  /r  ) /t  i  -  r  ■  / r  i  > 

In  the  limit  where  (  T  .  /  T  )>  I  a  dependence  ol  this  ratio 
on  temperature  w  ill  have  little  effect  on  k  .  but  in  the  genei al 
case  one  would  have  to  make  assumption-  concerning  it' 
variation  with  temperature  Here,  we  are  interi'ted  in  the 
temperature  dependence  of  the  ratio  ot  two  decay  mode'  and 
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it  seems  plausible  to  assume  a  temperature  dependence  of  the 
form 

(r,/r.)  =  (Tc/T)  -  \  (24) 

where  T  is  an  adjustable  parameter  to  be  found  by  compari¬ 
son  with  experimental  data  and  is  that  temperature  for 
which  *  =  0  5  kL  ,  i.e..  one-half  of  all  collisions  result  in 
stabilization.  The  choice  of  the  exponent  y  is  not  obvious.  As 
an  estimate  one  might  take>  as  equal  to  one-half  of  the  differ¬ 
ence  of  the  rotational  degrees  of  freedom  of  the  decay  prod¬ 
ucts  in  reactions  (21 )  and  (22).  Assuming  that  free  internal 
rotations  are  present  in  the  unstable  complex  but  not  in  the 
stable  molecule  and  that  the  latter  has  a  linear  structure  one 
obtains  >•  =  25  and  this  choice  was  adopted. 

In  order  to  learn  if  such  a  model  could  reproduce  the 
experimental  data,  three-body  rate  coefficients  were  calcu¬ 
lated  using  Eq  ( 8 )  and  taking  k  ( 8  )  =  k  as  defined  by  Eqs. 
(2.-)  and  t  24 ).  The  values  of  l'(  T)  at  T  =  300  K  were  ob¬ 
tained  by  interpolating  those  given  by  Bates,'  who  numen- 
caily  calculated  absolute  values  for  a  number  of  assumed 
potential  well  depths  of  the  product  ions  assuming  free  inter¬ 
nal  rotations  in  the  complex.  Adjustments  for  different  rota¬ 
tional  and  v  irrational  constants  were  made  using  the  inter¬ 
polation  taPle  given  by  Bates'  and  molecular  constants  from 
Herzberg  J  and  Albritton  el  al  :'  tfor  NO’  ) 

The  Langev in  (induced  dipole)  rate  was  used  for  kL  in 
Eq  23  taking  molecular  polarizabilities  from  the  table  giv¬ 
en  by  McDaniel  and  Mason  o'*  The  exponent  >•  was  taken  to 
have  a  value  of  2  5  for  all  reactions  such  that  T  remained  as 
the  only  parameter  that  depended  on  the  specific  reaction. 

The  results  of  the  model  calculations  are  compared  to 
the  experimental  data  in  Figs  1  to  4  In  all  cases  values  of  T, 
could  be  chosen  such  that  the  calculated  values  of  the  three- 
body  rate  coefficients  agree  quite  well  with  the  experimental 


r  l«  R*  .  *-tfuxnt  for  react;  *n  if  »'  a  tuiu  .1*  of  temper  at  sire 
.i'v  \  .  *  ;■ v-'inenfal  daU  Ref  !■  Da>h<d  line  mtxJel  cauula 


FIG.  2  Rate  coefficient  for  reaction  ( 1 7)  as  a  function  of  gas  temperature. 
Open  circles,  data  from  Ref  1 1  Squares  data  from  Ref  12  Dashed  line 
model  calculation 

data.  The  model  calculations  shown  in  the  figures  have  the 
functional  forms: 

k(T)  =  *(300)  (300/7V  [1  +  (T/Tc)2f]  (25) 

with  the  following  parameters: 

Reaction  (16),  *(300)  =  4.6x  10  ^cmVs,  Tc  =  1000K; 

reaction  (17),  *(300)  =  1.5x  lO'^cmVs,  Tc  =200  Kq 

reaction  ( 18),  *(300)  =  5.3 X  10  wcm7s,  Tc  =  150  Kj 

reaction  (19),  *(300)  =4.0x10  M’cm7s,  Tc  =  120  K. 

A  noticeable  deviation  from  a  simple  power  law  (“cur¬ 
vature"  in  log-log  plots  of  *  vs  7~)  is  present  in  the  case  of 
reaction  (2)  (see  Fig.  2).  The  low-temperature  data  by 
Rowe  et  al.':  and  Bohringer  and  Arnold"  seem  to  yield  a 
weaker  temperature  dependence  at  lower  temperatures,  but 


1  IG  3  Rate  x.'rfhcicni  for  nation  <  1 H  i  as  a  furu lion  of  temperature  Or- 
Jes  data  from  kef  Id  Data  p«»ini  I  .the  led  MHKR  from  Ref  P  Dashed 
line  rnixlel  t.ik ulalum 
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PIG  4  Rale  coefficient  for  reaction  *19)  as  a  function  of  temperature 
Open  circles  data  from  Ref  10  Points  labeled  JHB.  SAG.  and  H  +  V  are 
from  Refs  I.'.  IQ.  and  18.  respective!)  Dashed  line  mode!  calculations 

unfortunately  the  experimental  data  are  not  sufficiently  ac¬ 
curate  to  permit  the  statement  that  the  model  calculations 
"explain”  this  departure  from  a  power  law. 

It  is,  of  course,  not  surprising  that  one  can  reproduce  the 
experimental  data  at  one  temperature  by  choosing  a  particu¬ 
lar  value  of  Tt ,  but  it  seems  remarkable  that  one  choice  of 
this  parameter  yields  both  an  acceptable  fit  to  the  absolute 
value  of  the  rate  coefficient  and  the  temperature  depen¬ 
dence.  The  fraction  [ see  Eq  ( 23 )  ] 

/?  =  ( r */r  )/<  1  +  r/r  >,  (26) 

which  may  be  identified  with  the  usual  “complex  stabiliza¬ 
tion  efficiency,"  is  near  unity  for  reaction  ( 16)  (ft  =  0.95) 
at  300  K.  but  assumes  small  values  (ft  —  0.27,  0  15,  and 
0.06)  for  reactions  ( 17 ).  1 18),  and  (19),  respectively,  at  300 
K.  As  one  might  have  expected,  the  efficiency  of  complex 
stabilization  strongly  decreases  with  decreasing  stability  of 
the  product  ion.  Ai  very  low  temperatures  7?  approaches 
unity  for  all  reactions  and  the  three-body  association  rates 
should  then  vary  as  7'  .  This  model  prediction  still  awaits 

experimerral  resting 

IV.  SUMMARY  AND  CONCLUSIONS 

It  has  been  shown  that  the  mea-.iicd  temperature  de¬ 
pendence  of  several  three-body  as  >x  ation  reactions  can  be 
reproduced  by  modeling  the  complex  stabilization  as  pro¬ 
ceeding  through  a  super  complex.  Ft  is  clear  that  the  theory 
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as  presented  here  is  essentially  semiempirical  since  critical 
parameters  are  derived  from  experiment,  but  it  is  still  re¬ 
markable  that  a  minimum  of  plausible  assumptions  leads  to 
quite  acceptable  agreement  with  experimental  data. 

The  model  indicates  that  ion  association  rate  coeffi¬ 
cients  may  not  vary  in  the  form  of  simple  power  laws  over  a 
wide  range  of  temperatures  Specifically,  it  is  likely  to  be 
incorrect  to  extrapolate  experimental  data  to  temperatures 
much  below  the  range  covered  by  the  experiment 

Complex  stabilization  via  formation  of  a  supercomplex 
is  likely  to  be  a  poorer  approximation  in  cases  where  the 
third  body  is  an  atom  rather  than  a  molecule.  This  unfortu¬ 
nately  includes  the  case  of  helium,  the  buffergas  most  often 
used  in  flow  tube  measurements  of  three-body  association 
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APPENDIX  III 


The  Electron  Temperature  Dependence  of  the  Recombination 
of  Electrons  with  NO*  Ions 

J.  L  Dulaney.  M.  A.  Biondi,  and  R.  Johnsen 
Department  of  Physics  and  Astronomy 
University  of  Pittsburgh 
Pittsburgh.  Pennsylvania  15260 


ABSTRACT 

The  electron-temperature  dependence  of  the  dissociative  recombination  of  NO  ions  with 
electrons  has  been  re-determined  using  a  microwave  afterglow /mass  spectrometer  apparatus 
employing  microwave  heating  of  electrons  in  an  effort  to  resolve  the  discrepancy  between 
previous  microwave  afterglow  measurements  of  the  recombination  coefficient  a(NO  )  at  high 
electron  temperatures  made  by  Huang  et  al.  (a  ~  Te  3  )  and  the  trapped-ion  results  of  Walls 
and  Dunn  (a  ~  Te  *5).  It  is  found  that  the  electron  temperature  attained  by  microwave 
heating  in  the  experiments  of  Huang  eL  al.  is  lower  than  calculated  by  them  due  to  neglect  of 
the  effect  of  inelastic  collisions  of  electrons  with  minority  nitric  oxide  molecules  and  that  the 
assumption  of  spatially  uniform  electron  temperatures  is  not  warranted.  An  improved  data 
analysis  technique  is  described  which  allows  accurate  determination  of  the  recombination 
coefficient  at  elevated  electron  temperatures.  Using  this  analysis,  the  recombination  coefficient 
is  found  to  follow  the  relation. 


ot(NCT)  =  (4.2  ±  .2)xl0"7  (Te(K)/300)"<  ’5±oo**  cm3/s  , 


over  the  electron  temperature  range,  295K  <  Tt  <  4500K.  at  T.  =  Tn  =  295K.  The  afterglow 
data  now  support  the  trapped-ion  results  of  Walls  and  Dunn. 


To  be  submitted  to  Physical  Review  A 


DT /C 


